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We investigate the unpolarized virtual photon structure functions F2 {x , Q'^ , P'^) and 
^ I F2{x, Q'^,P'^) in perturbative QCD for the kinematical region < < Q^, where — P^) is 

' the mass squared of the probe (target) photon and A is the QCD scale parameter. Using the frame- 
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work of the operator product expansion supplemented by the renormalization group method, we 
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QJ [ derive the definite predictions for the moments of Pj(x,Q^,P^) up to the next-to-next-to- leading 
'fl ' 

^ I order (NNLO) (the order aog) and for the moments of F2{x,Q^ , P'^) up to the next-to-leading 
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^ ■ order (NLO) (the order aag). The NNLO corrections to the sum rule of F2{x, Q'^,P^) are negative 
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. . . , and found to be 7% ~ 10% of the sum of the LO and NLO contributions, when P^ = IGeV^ and 

q2 = 30 ~ lOOGeV^ or P^ = 3GeV^ and = lOOGeV^ and the number of active quark flavors 
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numerically to obtain the predictions for F2{x, Q^, P^) and F2{x, Q^, P^) as functions of x. 
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FIG. 1: Deep inelastic scattering on a virtual photon in the e+ e collider experiments. 

I. INTRODUCTION 

The experiments at the Large Hadron Collider (LHC) will be started shortly and it is 
much anticipated that signals for the new physics beyond the Standard Model (SM) will 
be discovered [1]. Once these signals are observed, they would be examined more closely 
in a proposed e+e~ collider machine called the International Linear CoUider (ILC) [2]. In 
analyzing these signals for the new physics, the knowledge of the SM, especially, of QCD will 
be more important than ever before. It is well known that, in e^e~ collision experiments, the 
cross section for the two-photon processes e^e~ — > e+e" -|-hadrons shown in Fig.l dominates 
at high energies over other processes such as the annihilation process e^e~ — > 7* ^ hadrons. 
Here we consider the two-photon processes in the double-tag events, where both of the 
outgoing e'^ and e" are detected. In particular, we investigate the case in which one of the 
virtual photon is very far off shell (large = — g^), while the other is close to the mass shell 
(small = —p^). This process can be viewed as a deep-inelastic electron-photon scattering 
where the target is a photon rather than a miclcon [3]. In the deep-inelastic scattering off 
a photon target, we can study the photon structure functions, which are the analogs of the 
nucleon structure functions. The photon structure functions are defined in the lowest order 
of the QED coupling constant a — e^/Air and, in this paper, they are of order a. 

The unpolarized (spin-averaged) photon structure functions F^{x,Q'^) and F2{x,Q'^) 
of the real photon (P^ = 0) were first studied in the parton-model (PM) [4] and then 
investigated in perturbativc QCD. A pioneering work was done by Witten [5] in which 
he derived the leading order (LO) QCD contributions to F2 and F^. A few years later 
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the next-to-leading order (NLO) corrections to were calculated [6]. These results were 
obtained in the framework based on the operator product expansion (OPE) [7] supplemented 
by the renormalization group (RG) method. The same results were rederived by the QCD 
improved PM powered by the parton evolution equations [8, 9]. Recently, the lowest six 
even-integer Mcllin moments of the photon-parton splitting functions were caluculated to 
the next-to-next-to-leading order (NNLO) and the parton distributions of real photon and 
the sturucture function were analyzed [10]. The same authors later gave the compact 
and accurate parameterization of the photon-parton splitting functions up to the NNLO in 
Rcf.[ll]. 

When polarized beams are used in e+e^ collision experiments, we can get information 
on the spin structure of the photon. The QCD analysis of the polarized structure function 
g1{x, Q^) for the real photon target was performed in the LO [12] and in the NLO [13, 14]. 
For more information on the theoretical and experimental investigation of both unpolarized 
and polarized photon structure, see Ref.[15]. 

A unique and interesting feature of the photon structure functions is that, in contrast 
with the nucleon case, the target mass squared —P^ is not fixed but can take various values 
and that the structure functions show different behaviors depending on the values of P^. 
The photon has two characters: The photon couples directly to quarks (pointlike nature) 
and also it behaves as vector bosons (hadronic nature) [16]. Thus the structure function 
F^(a:, Q^) of the real photon (P'^ — 0) may be decomposed as 

F^{X, Q^) = F^{X, g')|pointlike + F^{X, g')|hadronic ■ (1-1) 

The first term, a pointlike piece, can be calculated in principle in a perturbative method. 
On the other hand, the second term, a hadronic piece, can only be computed by some 
nonperturbative methods like lattice QCD, or estimated, for example, by the vector meson 
dominance model [16]. 

The moments of -^2^(2;, (5^)| pointlike and -^2^(2;, (5^)|hadronic for even n may be written, re- 
spectively, as 

2 dxx^-^F^{x,Q%,,^,i^, = a[-^an + bn + 0{as{Q'))] , (1.2) 

/ dxx"-'F2'^(x,Q')|hadronic = a/i„(a,(g')) , (1.3) 
Jo 
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where x is the Bjorken variable and agiQ"^) — g'^{Q'^)/A7r is the QCD running couphng 
constant. Since l/as{Q'^) behaves as ln((5^/A^) at large-Q^, where A is the QCD scale 
parameter, the first term an/agiQ"^) dominates over the 6„ term and also over the hadronic 
term hn{ois{Q^))- The definite prediction for the LO contributions a„ was given in Ref.[5]. 
Meanwhile, the NLO corrections 6„ were calculated only for n > 2 in Rcf. [6]. For n> 2, 
the hadronic moments hnioisiQ"^)) vanish in the large-Q^ limit and the 6„ terms give finite 
contributions. However, at n — 2, the hadronic energy-momentum tensor operator comes 
into play. Due to the conservation of this operator, 6„ shows a singularity at n = 2 and 
^n=2(as(Q^)) does not vanish at large-Q^. Actually, hnioisiQ"^)) also develops a singularity 
at n = 2 which cancels out the one of 6„, and hn{as{Q^)) and 6„ in combination give a 
finite but perturbatively uncalculable contribution at n = 2 [17]. The fact that a definite 
information on the NLO second moment is missing prevents us to fully predict the shape 
and magnitude of the structure function of F^{x, Q"^) up to the order 0{a). 

It was then pointed out in Ref.[17] that the situation changes significantly when we 
analyze the structure function of a virtual photon with much larger than the QCD 
parameter A^. More specifically, we consider the following kinematical region, 

A^ < p2 ^ (1.4) 

In this region, the hadronic component of the photon can also be dealt with perturbatively 
and thus a definite prediction of the whole structure function, its shape and magnitude, 
may become possible. In fact, the virtual photon structure function F2{x,Q'^, P'^) in the 
kinematical region (1.4) was calculated in the LO (the order a/ag) [18] and in the NLO (the 
order a) [17, 19], and the longitudinal structure function Fl{x, Q^, P^) in the LO (the order 
a) [17] without any unknown parameters. It is notable that the pathology of singularity, 
which appeared at n — 2 in the term 6„ of Eq.(1.2) for the real photon target, disappeared 
from the moments of P^(x, Q^, P^). The parton contents of the virtual photon for the case 
(1.4) were studied in Refs. [20-22] . 

In the same kinematical region (1.4), the polarized virtual structure function gj{x, Q^, P^) 
was investigated up to the NLO in QCD in Ref.[23] and in the second paper of [14]. More- 
over, the polarized parton distributions inside the virtual photon were analyzed in various 
factorization schemes [24]. Quite recently the first moment of gJ{x,Q'^, P'^) was calculated 
up to the NNLO [25]. 
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FIG. 2: Forward scattering of a virtual photon with momentum q and another virtual photon with 
momentum p. The Lorentz indices are denoted by v, p, r. 

In this paper we investigate the unpolarized virtual photon structure functions 
(5^, P^) and F2{x, Q'^, P^) in the kinematical region (1.4) in QCD. Here we neglect all 
the power corrections of the form {P"^ /Q'^Y = 1; 2, • • • ) which may arise from target mass 
and higher-twist effects. We present definite predictions for F^(a;, Q^, -P^) up to the NNLO 
(the order aag) and for F2{x,Q^,P'^) up to the NLO (the order aas). The recent calcu- 
lations of the three-loop anomalous dimensions for the quark and gluon operators [26, 27] 
and of the three-loop photon-quark and photon-gluon splitting functions [11] have paved 
the way for this investigation. Using the framework of the OPE supplemented by the RG 
method, we give, in the next section, an expression for the moments of F^ix, Q^, P^) up to 
the NNLO corrections. In Sec. Ill we enumerate all the necessary QCD parameters to evalu- 
ate the NNLO corrections. In Sec. IV the second moment of -^^(a;, Q^, P^) will be evaluated 
up to the NNLO. The numerical analysis of F2{x, Q^, P^) as a function of x will be given 
in Sec.V. In Sec. VI the longitudinal virtual photon structure function Fl{x, Q^, P^) will be 
analyzed up to the NLO. The final section is devoted to the conclusions. 

II. THEORETICAL FRAMEWORK BASED ON THE OPE AND THE NNLO 
CORRECTIONS TO F^{x,Q^,P^) 

In this article we analyze the virtual photon structure functions F2{x,Q^, P"^) and 
Fl{x,Q'^,P'^) using the theoretical framework based on the OPE and RG method. Un- 
less otherwise stated, we will follow the notation of Ref.[6]. Let us consider the forward 
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virtual photon scattering amplitude for 7(g) + j{p) — > 7(g) + 7(p) illustrated in Fig.2, 

T^uAp, <l)-i j d^xd'yd'ze'^-^e'P<^-'\0\T{J^{x)MO)J,{y)Mzm , (2.1) 

where is the electromagnetic current. Its absorptive part is related to the structure tensor 
Wfj^uprip, q) for the target photon with mass squared — —P^ probed by the photon with 

W^^prip, q) = -ImT^^^ip, q) . (2.2) 

71 

Taking a spin average for the target photon, we get 

w;i^{p,q) = ^E<A)(p)^.-pr(p,g)6[,)(p) 

A 



—g'^w^uAp^q) 



^^J c^^xe^^^(7(p)|JM(x)J.(0)|7(p))spinav (2.3) 

Now WjJ^,{p,q) is expressed in terms of two independent structure functions F2{x,Q'^ , P"^) 
and F^{x,Q^,P^): 

W^p^q) = {5..-^}^F2(x,Q^P^) 

.{-..14^- (^^1>'^'^^^' ^'-'^ 

where x = Q'^/2p ■ q. 

Applying OPE for the product of two electromagnetic currents at short distance we get 

n=even 

+ \-9n\9uaq^ + Qnxq^qa + 9uaq,j,q\ - 9nyqxqo\ 

E (^)?W-?M„-2E^2,nO^ 



n=2 
n=even 



+•••, (2.5) 

where and C2,n coefficient functions which contribute to the structure functions 

F2 and F^, respectively, and Of^'"^"' and q^'^i^^''''^^-^ spin-n twist-2 operators (hereafter 
we often refer to O^i "''" as Of). The sum on i runs over the possible twist-2 operators 
and • • • represents other terms with irrelevant coefficient functions and operators. In fact, 
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the relevant Of are singlet quark (■0), gluon (G), nonsinglet quark (A^-S") and photon (7) 
operators as follows: 

QW-Mn ^ i^-^^^it^^D^"^ ■ ■ ■ Di'^'^ lip - trace terms , (2.6a) 

O!^'"'"" = - i^'-'^Gj^'D^'^ ■ ■ ■ DMn-i(;;aMn} _ tj-j^ce terms , (2.6b) 

O^i^"'^" = i"-iV^7^^iL)'^2 . . . D^"^{Ql^ - {e^)l)il: - trace terms , (2.6c) 

QMi-Mn ^ ^ i^^-^Fj^W^^ ■ ■ ■ a^"-iF"''"> - trace terms , (2.6d) 

where { } means complete symmetrization over the Lorentz indices //i • • • //n and D^^ denotes 
covariant derivative. In quark operators and O'^g given in Eqs.(2.6a) and (2.6c), 1 is an 
Uf X Uf unit matrix, Q^f^ is the square of the Uf x Uf quark-charge matrix, with Uf being 
the number of active quark (i.e., the massless quark) flavors, and (e^) = (X^r^ ^D/^f is the 
average charge squared where Cj is the electromagnetic charge of the active quark with flavor 
i in the unit of proton charge. It is noted that wc have a relation Tr{Q^i^ — (e^)l) = 0. The 
essential feature in the analysis of the photon structure functions, in contrast to the case 
of the nucleon counterparts, is the appearance of photon operators O" in addition to the 
familiar hadronic operators O^, Oq and O^^^^ [5]. 

The spin-averaged matrix elements of these operators sandwiched by the photon state 
with momentum p are expressed as 

(7(p)|0^'^"■•'^17(p))spin av = A^f,' , P') {p^^ ■ ■ ■ p^^ - tracc terms} (2.7) 

with i — ip,G, NS, 7, and /i is the renormalization point. Then the moment sum rules for 
and F2 are given as follows [7]: 

' dxx^-'F^(x, Q\ P') = Yl <^2,n(QV/^', 9(^^% «) Aii^i', P') , (2.8a) 

i=ilj,G,NS,'y 

\xx--'Fl{x,Q\P^) = Yl Q,n(QV/^',^(/^'),«)^U/^',^'), (2.8b) 

i=ip,G,NS,'r 

with gifi"^) being effective running QCD coupling constant at yU^. Recall that in this article 
the photon structure functions are defined to be of order a. Since the coefficient functions 
„ and C2,n ^(o^); it is sufficient to evaluate at 0{1). Thus we have 

Al{i,\P') = l. (2.9) 
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On the other hand, the matrix elements {i — t/j,G, NS) for the hadronic operators start 
at 0{a). For —p^ — P'^ ^ A^, we can calculate {i — ip^G^NS) perturbatively in each 
power of g'^. When /x^ is chosen at P^, they are expressed in the form as 



Ai{^^', P')\,.=p2 = ^A\i9{P')) , for i = V', G, NS . 



(2.10) 



Let us first analyze the structure function (x, Q^, P^). We will evaluate its moment sum 
rule up to the NNLO. The Q^-dependence of the coefficient functions C2^^{Q'^ / jj? ,g{iJ?),a) 
in (2.8a) is governed by the RG equation. Putting //^ = —p^ — P^, its solution is given by 



C^,jg7P^^(P^),a)= Texp 



7n(fl',a) 



-9{Q') M 



CiJirg{Q'),a) , (2.11) 



with i, J = '0, G, NS and 7 . Here (3{g) is the beta function and ^n{g^-i Oi) is the anomalous 
dimension matrix. To the lowest order in a, this matrix has the following form: 



7n («/,«) 



In 



K„{g, a) 







(2.12) 



where 7„(5'^) is the usual 3x3 anomalous dimension matrix in the hadronic sector 



inig) = 



, ^^^s{g)) 

and Kn{g, a) is the three-component row vector 



(2.13) 



Kn{g,a) = (^K'^{g,a) , K'^{g,a) , K'^gig^^^)) > 



(2.14) 



which represents the mixing between photon operator and remaining three hadronic opera- 
tors. Then the evolution factor in (2.11) is expressed in the form as [6], 



Texp 



9(Q2) Pig) 



M„ 










(2.15) 



where 



M^{QyP^,g{P^)) = Texp 
X^{QyP',g{P'),a) = 



^1 



-gm P(g) 



dg r,, X ^exp 

g(Q2) P\g) 



g{Q'^) P\g ) 



(2.16) 
(2.17) 
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Thus using (2.9)-(2.11) and (2.15)-(2.17), we get 



a 



dxx-'F^{x,Q',P') = ^A^{-g{P'))-M^{Q'/P^-giP')).C2,n{h9{Q')) 



with 



and 



+x„(gvp^^(p^),a) .C2,„(i,^(g^)) 

^n{9)-{At{g)J'i{g)Al\9)). 



C2,n(l,^) 



(2.18) 



(2.19) 



(2.20) 



In order to evaluate Mn{Q'^/P'^,g{P'^)) in (2.16) up to the NNLO, we first expand %{g) 
in powers of g'^ up to the three-loop level as 



„2 4 6 

167r2 ^" ^ (167r2)2 '^^ ^ (167r2)3 



(2.21) 



Then putting gi = g{P^) and g2 = giQ"^), wc find that Mn{Q'^/P'^,g{P'^)) is expanded as 

r 

M^{QyP^g{P'))^Texp / dg- 

-J 92 



Ma) 



— exp 



9i 



92 



Pig) 



91 



dg exp 



92 



91 



+ / dgexp 

'92 



r9i r9] 
+ / dgaBxp 1 

Jg2 J 9a 



^ln\gb) 
X— -— exp 

PKQb) 



9b 



dg" 



92 



M 
jlil) 

m 

li'\9a) 

Piga) 

P{9"') 



exp 



exp 



92 



9 ^'^^^(n"\ 



92 



Pig") 



9a 



dgb exp 



92 



' Pig") 



9b 



+ 



(2.22) 



To evaluate the integrals, we make a full use of the projection operators obtained from the 
one-loop anomalous dimension matrix 7^°^ in (2.21) [6]: 



9(0) 

In 



i=+,-,NS 



(2.23) 
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where A" (i — +, —, NS) and P" are eigenvalues of 7^*^^ and the corresponding projection 
operators, respectively. The explicit forms of A" and are given in Appendix A. Expanding 
(3{g) in powers of up to the three-loop level as 

(2.24) 

we perform integration in (2.22). The final form of M^iQ'^ / P'^ , g{P'^)) up to the NNLO is 
given in (A7) in Appendix A. 

Similarly, expanding Kn{g, a) in powers of g^ up to the three-loop level as 



g'^ 

^^^^ ^ ~16^^' ~ Jl6^^' ~ (16^^^^ + 



Kn{g,a 



\ -*- ^ r> 



16n^"'' (167r2)2 

2 / d2 d2 



o2„4 



^(1) _ e g (,) 
" (167r2)3 " 



(2.25) 



we can evaluate Xn{Q^/P^,g{P^), a) in (2.17) up to the NNLO. The result is given in (A8) 
in Appendix A. 

Finally, expansions are made for the photon matrix elements of hadronic operators 
AnigiP"^)) in (2.19) as well as the coefficient functions C2,„(l, ^(Q^)) in (2.20) and 

C2 „(1, ^(Q^), a) in (2.18) up to the two-loop level as follows: 



AMP')) - A^:^ + ^^i'^ + 



12^^ _ ^(0) gHQ^^ii) 9\Q'^) ^(2) 



C2,n(l,^(g^)) = cfu 



cU-^r9{Q%a) 



e-'fiQ') 



(2.26) 
(2.27) 

(2.28) 



167r2^^'" ' (167r2)2 ^^'^ 
Then putting (A7), (A8), (2.26), (2.27) and (2.28) into (2.18), we obtain the expression for 
the moment sum rule of F^{x, Q^, P^) up to the NNLO (aag) corrections as follows: 



^0 



dxx''-^F^{x,Q'',P'' 



a 1 



47r 



47r2/3o (q;,(Q2)Z. 

i 



1 - 



as{P^ 



»s{P'), 



1 - 



asiP' 
as{P')J 



d"-l 



i 

-Q 



as{P') 
1 



as{P' 



+0{al 



with i = +,-,NS, 



(2.29) 
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where (i? = ■ The coefficients O^, A^, Bf, C^, V^, and are given by 



^ = ^"C-S^ ' (2-30) 



-rK^^Prcfl^ - 2/3oAW i^-C^S , (2.31) 

wn ^ ir^(O) > ^3 ^(0) 

r = 2/3o(<^ + 1« • eg) , (2.33) 

pnx,(2) pn 1 



1^(0) -^fc 7« -r:,- In (0) i 

" 2^ (\n _ \n ,cyn \(^n _ \n , ARA ^ 



- (A« - A- + 2/3o) (A- - A" + A(3o) ^''^ 1 - 



o pn^(l) pre -1 

, pn^{0) Pi 1^(1) /" ^(0) ^ 



-2/3oAWi^«c(°)^cir - 2/3oAf i^"Cg , (2.34) 



11 



cn ^ 1^(0) p«/-^(l) A 1 ~ 1^(0) -Pj"^" ^(1) 1 



/Q2 



_ 1^(0) -P"7n -P" (0) /3i 1 - (i" P-'^i^'Pp „(0)/3i 

" ^Af-A7 + 2/3o ''"/3o < " ^A"-Af + 2/3o ''"/3o 

pn - (1) pn - (1) pn i 

" f^(A?-A« + 2/?o)(A«-Af + 2/?o) 

- /i-^l) p«r(°) — + ft-^l) V -^"7n ^Pj" ^(0) J_ 

pnr>{l) pn n 



(2.35) 



-rn _ 1^(0) pn^(2) 1 1^(0) pn,^(l) A 



pn^X'-) pr. 

^A"-A'^ + 2/3o + 



/3( /32 1 A rfj' 

2 




I ix-(O) pn>-((0) 



pn^(2) pn -. 

I 1^(0) M 7n (0) 1 

" VAf-A^" + 4/3o "'"l + df 

, i>-(0) \^ ^ i ^" ^^i ^" " k ^(0) 
n Mn _ xn , OR^"^' 



J 



-Tj 7n -T-/ ^/n -r,, 



-(A^--A- + 2;3o) ^■nAf-A^" + 2/3o A" - A^ + 4/3o7 1 + ci^ 

(2.36) 

= 2/3o(c5 ^ + A« • C(^) + li^) • C^S) , (2.37) 

with ij,k = +, -,A^5. The LO term was obtained by Witten [5]. The NLO (a) 
corrections A^, B" and C" without terms with A^^ were first derived by Bardeen and 
Buras [6] for the case of the real photon target (i.e. = 0). Later authors in Ref.[17] 
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analyzed the NLO (a) corrections for the case of the virtual photon target (P^ ^ A^) and 
the terms with An"* were added to Af and C". The coefficients Vf, and Q'^ are the 

NNLO {aas) corrections and new. 

For n = 2, one of the eigenvalues, A"^^, in Eq.(2.23) vanishes and we have dPz""^ = 0. 
This is due to the fact that the corresponding operator is the hadronic energy-momentum 
tensor and is, therefore, conserved with a null anomalous dimension [6]. The coefficients 
A":r'^ and have terms which are proportional to and thus diverge. However, we 

see from (2.29) that these coefficients are multiplied by a factor 1 — (as{Q'^)/as{P^) ) 
which vanishes. In the end, the coefficients A'^"^ and multiplied by this factor remain 
finite [17]. 

III. PARAMETERS IN THE MS SCHEME 

All the quantities necessary to evaluate the NNLO (aag) corrections to the moments 
of F^{x, Q^, P^) have been calculated and most of them are presented in the literature, 
except for the two-loop photon matrix elements of hadronic operators An^^ , An^'^ and 
^(j2)NS ^ Also for the three-loop anomalous dimensions K^^'"^, Kq^'"' and K^^^g^ , we only have 
approximate expressions in the form of photon-quark and photon-gluon splitting functions. 
In the following we will enumerate all these necessary parameters. The expressions are the 
ones calculated in the modified minimal subtraction (MS) scheme [28]. 

A. Quark-charge factors and f3 function parameters 

The following quark-charge factors are often used below: 

Uf Uf 

= (e') = E e'/^/ ' ^NS^l, S,^ 3nf{e') = 3 . (3.1) 

i=l i=l 

The P function parameters /3o, Pi and P2 [29] are given by 

Po = - \nf , (3.2) 

A = ^-^C\ - ^^CAUf - 2CFnf , (3.3) 

2857^3 1415^2 205^ ^ 79^ 2 ^2 2 

= - -^Clrif - —CACFUf + -Cau) + Cln^ + -C^n) , (3.4) 

with Ca = 3 and Cf = | in QCD. 
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B. Coefficient functions 



As shown in (2.27) and (2.28), we need the hadronic coefficient functions C2_„(l, ^(Q^)) 
with i — t/j,G and NS, and the photon coefficient function C2^^{l,g{Q'^),a) up to the 
two-loop level. At the tree-level, we have 

The one-loop coefficient functions were calculated in the MS scheme in Refs.[28] and [30]. 
The MS results are written in the form as 

Ct,n^ = S^B^, C^^^^ = 8^Bq , C^f = ^NsBj^g , Cl^l'' = S^B^ , (3.6) 

where = ^ns ^^'^ obtained, for example, from the MS-scheme results for = 

B'^g and Bq given in Eqs.(4.10) and (4.11) of Ref.[6] by discarding the terms proportional 
to ln(47r — 7^;). B^ is related to Bq by 'B^ = {2/nf)'B^ . 

The two-loop coefficient functions corresponding to the hadronic operators were calcu- 
lated in the MS scheme in Ref.[31, 32]. They were expressed in fractional momentum space 
as functions x. The results in Mellin space as functions of n are found, for example, in 
Ref. [33]: 

CtS^ = SAcfl^^'^'in) + 4J-(n) + c2'-(n)}, (3.7) 
= S^cfhn), (3.8) 



= ^^.{4r^^(n) + c2'-(n)}. (3.9) 



I 

where c^J'^°'(n), c^J'"'''(n), c^J''''(n) and 4^(n) are given in Eqs.(197), (198), (201) and 
(202) in Appendix B of Ref. [33], respectively, with N being replaced by n. The two- loop 
photon coefficient function 6*2^^^ is expressed as 

<^ = sA^in) , (3.10) 
where cf}^{n) is obtained from C2^p(n) in (3.8) by replacing Ca — >■ and ^ — >• 1 [10]. 

C. Anomalous dimensions 

The one-loop anomalous dimensions for the hadronic sector were calculated long time 
ago [34, 35]. The expressions of 7^°^'" = Ins^^ T^^g'"' 7g|'" 7gg" given, for example, 
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in Eqs.(4.1), (4.2), (4.3) and (4.4) of Ref.[6], respectively, with / being replaced by Uf. As 
for the one-loop anomalous dimension row vector Kn^ — (K^^'^ , Kq^'"' , K^^^g^), we have 
Xg)'" = 0, and irf and K^^]^'' are given, respectively, in Eqs.(4.5) and (4.6) of Ref. [6] 
with / being replaced hj nj again. 

The two-loop anomalous dimensions for the hadronic sector were calculated in Ref. [30] 
and recalculated using a different method and a different gauge in Ref. [36] . The results by 
the two groups agreed with each other except in the part of 7^^" proportional to Cq, but 
this discrepancy was solved later [37]. They are given by 



(l),n 

Tns 


= 27(^,)+(n) , 


(3.11) 


(l),n 


= 2(7(^.)+(n)+7«(^)) , 


(3.12) 


(l),n 
Ii/jG 


= 27SHr^) , 


(3.13) 


(l),n 


= 27W(n) , 


(3.14) 


(l),n 

7gg 


- 27S)(r^) , 


(3.15) 



where 'yns'^{n) is given in Eq.(3.5) of Ref.[26], and ^^\n), 7qg (n), 7gq (n) and 7gg (n) are 
given, respectively, in Eqs.(3.6), (3.7), (3.8) and (3.9) of Ref. [27], with N being replaced 
by n. The factor of 2 in (3.11) - (3.15) appears since, in Refs.[26] and [27], the anomalous 
dimension 7 of the renormalized operator O is defined as dO/dhiii^ = — 7O instead of 
(iO/ciln/x= -7O. 

The two- loop anomalous dimensions K^^'^, K^jPg^ and Kq^'^ can be obtained from 7^cr'" 
and 7gg" by replacing color factors with relevant charge factors [6]. Moreover we need an 
additional procedure for Kq^'^. They are given by 

X^'^'" = -3nf{e'')CFD^G{n) , (3.16) 
= -Srifiie') - {e^f)CFD^G{n) , (3.17) 
= -2>nf{e')CF{DGG{n) - 8) , (3.18) 

where D^ci^) and Dgg{^) a-i'c obtained from 7|/(^'" and 7^^^"' respectively, by replacing 
Ca and CpUf 2. The number 8 in (3.18) is due to the gluon self-energy contribution 
to 7gg", which should be dropped for X^^'" [38, 39]. 

The three-loop anomalous dimensions for the hadronic sector have been calculated re- 
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cently in Refs.[26] and [27]. They are expressed as 

tS" = 27(r(n) , (3.19) 

if/ = 2(7(rW+7gH^)), (3.20) 

7^'" = 27^)^ , (3.21) 

7^'" = 27£)(n) , (3.22) 

7^" = 27g)(n) , (3.23) 

where 7^^'''(n) is given in Eq.(3.7) of Ref.[26], and 7ps^(n), ^'^{n), 7gq (n) and 7gg^(n) are 
given, respectively, in Eqs.(3.10), (3.11), (3.12) and (3.13) of Ref.[27], with N being replaced 
by n. 

Concerning the three-loop anomalous dimensions Kf'''\ A'^^'" and xj^^'^ , the exact ex- 
pressions have not been in literature yet. In fact, the lowest six even-integer Mellin moments, 
n = 2,---,12, of these anomalous dimensions were calculated and given in Ref.[10]. Quite 
recently, the authors of Ref.[10] have presented compact parameterizations of the three- loop 
photon-non-singlet quark and photon-gluon splitting functions, Pnsj{x) and Pg^\x), instead 
of providing the exact analytic results [11]. It is remarked there that their parameterizations 
deviate from the lengthy full expressions by about 0.1% or less. They also gave in Ref.[ll] the 
analytic expression of the three- loop photon-pure-singlet quark splitting function pIs\{x). It 
is true that we can infer the analytic expressions for some parts of K^'^, K^g"' and Kq^'" 
from the known three-loop results of 7^^"^ and ^Qff'. For instance, the expressions of K^'" 
and K^^s^ which have the color factor Cp are obtained from 7^^'" by taking the terms which 
are proportional to the color factor rifCp. Also the terms of Kq^'^ which have the color 
factors UfCp and Cp are related to the ones of 7gg" with the color factors njCp and UfCp, 
respectively. But at present we do not have the exact analytic expressions of Kf'"", K^'g" 
and ir^^^'" whole. 

Under these circumstances we are reconciled to use of approximate expressions for K^'^, 
^Ns"^ and K^^'". They are obtained by taking the Mellin moments of the parameterizations 
for P^s\{x) and Pgj\x), and of the exact result for P^s\{x) which are presented in Ref.[ll]. 
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TABLE I: Numerical values of E^^gj{n), E^^y^"^{n), EQ^{n) and E'^^^^^ [n) for the lowest six even- 
integer values of n. The values for E^g^{n) {EQ^{n)) are found in Eq.(3.1) (Eq.(3.3)) or obtained by 
evaluating Eqs.(A.l)-(A.6) (Eqs.(A.7)-(A.12)) of Ref.[10]. The values of £;n^^'"°''(n) sxid E^^'"'' {n) 
are obtained from the expressions given in (B2) and (B3) in Appendix B, respectively. 













EG^{n) 




n = 


2 


- 86.9753+1.47051 




- 86.9844+1.47104 




31.4197 +5.15775 n/ 


31.4155 +5.15803 n/ 


n = 


4 


-102.831 +1.47737 




-102.848 +1.47787 


^/ 


23.9427 +1.10886 nf 


23.9419 +1.10888 n/ 


n = 


6 


-109.278 +1.65653 


n/ 


-109.299 +1.65699 


nf 


15.6517 +0.695953n/ 


15.6507 +0.695944n/ 


n = 


8 


-111.167 +1.69550 


nf 


-111.192 +1.69592 


nf 


10.9661 +0.498196n/ 


10.9651 +0.498178n/ 


n = 


10 


-111.035 +1.67061 


Uf 


-111.062 +1.67099 


nf 


8.16031+0.379060n/ 


8.15953+0.379038n/ 


n = 


12 


-109.943 +1.61908 


nf 


-109.972 +1.61943 


nf 


6.34829+0.300274n/ 


6.34777+0.300250n/ 



Then we have 



Xjj)'" ^ xS"pprox ^ -3n/((e^) - (e^)^) 2E^^^'^^-{n) , (3.24) 
i^f'" ^ ^^a^prox ^ -3n,(e2)2{E— (n) + E,,,{n)} , (3.25) 
^g),n ^ ^g),n^^^^ ^ -3n,(e2)2£;— (n) , (3.26) 

where the exphcit expressions of E^^^^°^{n), Eps-yin), and i?^'^^'^™(n) are given in Appendix 
B. Again the appearance of the factor of 2 in (3.24) - (3.26) is due to the difference in 
definition of the anomalous dimensions. As mentioned earher, the lowest six even-integer 
MeUin moments, n — 2, - ■ ■ , 12, of -f^^^s", K^^''^ and X^^'" were given in Ref.[10]. When we 
write K^^^g^ and Kq^'^ as 

Kp^n ^ _srif{{e') _ (e2)2) 2E^,,{n) , (3.27) 
Kg^'- = -3nf{e^)2EG^{n) , (3.28) 

then we get the exact results of Ensyin) and Ecyin) for even n = 2, ■ ■ ■ ,12. We give in 
Tablel the results of E^s-yin), E^^^^°^{n), EQ^{n) and £'Q^^™^(n) in numerical form for the 
lowest six even-integer values of n. We sec the deviations of E^i^'^°^[n) from E^^{n) and 
Eq^^^'^I'ti) from EQ^{n) are both far less than 0.1% for these values of n. 
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D. Photon matrix elements 



The two-loop operator matrix elements have been calculated up to the finite terms by 
Matiounine, Smith and van Neerven (MSvN) [40]. Using their results and changing color- 
group factors, we obtain the photon matrix elements of hadronic operators up to the two- loop 
level. 

First wc clear up a subtle issue which appears in the calculation of the photon matrix 
elements of the hadronic operators. The one-loop gluon coefficient function in (3.6) was 
calculated by two groups, BBDM and FRS (we have taken initials of the authors of Refs.[28] 
and [30] , respectively) . Both groups evaluated one-loop diagrams contributing to the forward 
virtual photon-gluon scattering as well as those contributing to the matrix element of the 
quark operator between gluon states, and they took a difference between the two to obtain 
B'^q. But actually BBDM calculated the gluon spin averaged contributions, i.e., multiplying 
QpT and contracting pairs of Lorentz indices p and r, whereas FRS picked up the parts which 
are proportional to QpT. Thus the BBDM results on the contributions to the forward virtual 
photon-gluon scattering and the gluon matrix element of quark operator are different from 
those by FRS, but the difference between the two contributions, i.e., S^, is the same as it 
should be. 

We have defined the photon structure functions F2 and in (2.3) and (2.4), taking a 
spin average of the target photon for the structure tensor W^ypr{p,q). We, therefore, adopt 
the BBDM result rather than that of FRS and convert it to the photon case. Then for the 
photon matrix elements of the hadronic operators at one-loop level we get , 

li^)^ = 3n^(e^)i/«(n) , 

IW^ = , (3.29) 

= 3n,((e^)-(e^)^)//«(n), 

where 



H?{n) 



11 4 4 /I 2 2 , ^ , , 

"~ + ~ ~ 7 — + ( A o\2 + + — ''^i(^) 

n (n + 1)"^ (n-|-2)"^ \n n + 1 n + 2 



(3.30) 



with Si{n) — Yl]=i 7 ■ Actually, Hq^\n) is related to the BBDM result on the one-loop gluon 
matrix element of quark operator A^^^ given in Eq.(6.2) of Ref.[28] as A^^q' — ^^H^^\n). 

MSvN have presented in Appendix A of Ref.[40] full expressions for the two- loop cor- 
rected operator matrix elements which are unrenormalized and include external self-energy 
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corrections. The expressions are given in parton momentum fraction space, i.e., in 2;-space. 
Taking the moments, the unrenormahzed matrix elements of the (flavor-singlet) quark op- 
erators between gluon states are written in the form as (see Eq.(2.18) of Ref.[40]), 



A. 



where 



^0 



dzz 



n—l Ak 
^9 



-p2 1 



) , k = PHYS , EOM and NGI 



(3.31) 



(3.32) 



and the expressions of A^hys(^^ i)^ a^om^z^_^^1) and IjGi(^, ^, 1) are given in 
Eqs.(A7), (A8) and (A9) of Ref.[40], respectively. Refer to Ref. [40] for the explanation of 



the "PHYS", "EOM" and "NGI" parts. The tensors T, 



(i) 



1, 2, 3) are given by (see 
Eqs.(2.19)-(2.21) of Ref.[40] and note that we have changed the Lorentz indices of gluon 
fields from /ii/ to pr), 

{A-pr, (3.33) 
{A-pr, (3.34) 



rp{2) ^ 

pr 



9pT - 
PpPr 



PpAr + ApPr ^ ApArP^ 

A-p {A-pyi 

PpAr + ApPr ApAy^ 



P^ 



pr 



A-p 

+ 



PpA, + App, , ApA^ 



{A-pf. 
A-pr, 



(3.35) 



2A-P '{A-pf. 

where is a lightlike vector (A^ — 0). The renormalization of Aqg^pT-{n, ^) proceeds 
as follows: First the coupling constant and gauge constant renormalizations are performed. 
Then the remaining ultraviolet divergences are removed by multiplication of the operator 
renormalization constants. We get the finite expression at n"^ — —p^ as 



47r/ 
An) 



"SWTS' 



pr 



(3.36) 



The expressions of aqg{n) and b, 



qg ' 



n] 



1,2) are given in Appendix C, while ^^^(n) is 



made up of the terms proportional to Ca-y ^''^^ is, therefore, irrelevant to the photon matrix 
element of the quark operator. Now multiplying g'^'^ and contracting pairs of indices p and 
T, we get 



1 



,pr. 



1 



-A 



2^' {A-py^"'''^ 



+ 



as 
Atx 
/a 



v47r 
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''qg 



^ (2)/ ■ 



}. (3.37) 



We can see from the expressions of aqg{n) and hqg{n) in (C2) and (C3), respectively, that the 
FRS result for the one-loop gluon matrix element of quark operator corresponds to Oqg {n) , 
while the BBDM result corresponds to the combination ja^g (n) — l^^g Indeed we find 

that H^^\n) in (3.30) is written as ^Hl''\n) = |aS(n) - ^^Sw}- 

The two-loop photon matrix elements of the quark operators are derived from the com- 
bination |aqg^(n) — ^b'qg{n) — in (3.37) with the following replacements: Ca — 0, 

and Cp x charge factor . The terms proportional to (^-^^ in Og5^(^) 

(2) 

and bqg (n) come from the external gluon self-energy corrections and should be discarded for 
the photon case. Thus we obtain 

= 3n,(e^)iff (n) , 

= 3nf[{e')-{er)H('Hn), (3.38) 



where 



-16S2,i(n)-48C3) 

_l_C ( \2 f ^ 



+-5i(n)( 



'6 8 16 _ 16 40 _ 32 

n + 1 ^ n + 2 ~ ^ ^ (n + 1)2 ~ (n + 2) 

4 80 56 16 48 64 

+ 



n n + 1 n + 2 n^ {n + iy {n + 2) 



2 



32 176 128 



n^ {n + lf {n + 2f 

/6 _8 16 ^0 32 >. 

^^^'^An^ n+l n^^ {n + lf {n + 2f) 

38 70 56 56 198 144 22 

^n n + l^ n + 2^ n^ {n + lf^ {n + 2y n^ 

40 128 20 88 1 

+ ZI + ■ (3-39) 



(n + l)3 (n + 2)3 n^ (n + 1)^ 

Similarly the renormalized matrix elements of the gluon operators between gluon states 
at /x^ = — p2 are written as (the unrenormalized version is given in Eq.(2.33) of Ref.[40]), 

2 
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The one- loop results a^g (n), bgg(n) and a^gl{n) are all proportional to the color factor C/ 



and thus they are irrelevant to the photon matrix elements. Also the two- loop result a^^' 
is made up of the terms proportional to C\ or and is irrelevant. The expressions 

(2) (2) 

of agg{n) and hgg{n) arc given by (C6) and (C7), respectively, in Appendix C. Then, we 
take the combination ^afg{n) — \hfg{n)^ and make replacements, Ca 0, {~2^ ~^ 
and ^ Cp x charge f actor j . Furthermore, we realize that the last two terms in 

parentheses of (C6) have also resulted from the external gluon self-energy corrections and 
are thus irrelevant for the photon case. In the end we obtain for the photon matrix elements 
of the gluon operators, 

If^ = 3n;(e^)4^)(n), (3.41) 

where 

■'\V ^ VV3(n-l) n n + 1 3(n + 2) {n + l)^J 

^ , , / 32 32 32 32 32 8 

+Si{n) -— + —— + - — — + ^ + 



+^-2(n)(. 



9(n-l) n n + 1 9(n + 2) v? (n + 1)2 

64 32 48 

~3(n + 2)2 ^ ^ ^ (n + 1)2 

32 32 32 32 \ 
h 



.3(n-l) n n + 1 3(n + 2)/ 
872 80 16 856 40 104 64 44 

^27(n- 1) ~ nT ^ n+1 ^ 27(n + 2) ~ ^ ^ (n + l)^ ^ 9(n + 2)2 ^ ^ 

+ 12 ^\ (3 42) 

(n + l)3 3(n + 2)3 n^ (n + 1)4 / ' ^ 

With all these necessary parameters at hand, we are now ready to analyze the moments 
of F2^(x, p2^ ^j^g NNLO. First we evaluate the coefficients C^, A^, Bf, C", Vf, S^, 
J^P and Q'^ with i — +, —,NS, the expressions of which are given in Eqs.(2.30)-(2.37), for 
n = 2,4,-- - ,12 in the cases of Uf — 3 and n/ = 4. The results are hsted in Tables II (for 
n/ = 3) and III (for nj = 4). In Table 1 of Ref.[17], the numerical values of the seven NLO 
coefficients, A'^, Bf, and C" with i = +, — , NS for n = 2, 4, ■ ■ ■ , 20 in the case of n/ = 4 
were already given. Our results on A^, and C" in Table III are consistent with those in 
Ref.[17] except for the values of A^ and A^. The discrepancy in the values of A^ and A^ 
arises from a term —8 in the parentheses of Eq.(3.18). See the discussion below Eq.(3.18). 
The numerical calculation of the NNLO coefficients and V'fifg for n = 2, 4, • • • ,12 in 

Tables II and III was performed by using the "exact" values of the three-loop anomalous 
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dimensions, kP^"", Kf'"^ and K^^^'", for n = 2, • • • , 12 given in Ref.[10] (at upper levels) 
and also by using the approximate expressions, K^Ns\pwo^^ ^i^approx and i^c approx defined 
in Eqs.(3.24)-(3.26) (in parentheses at lower levels). The coefficients .4." and cannot be 
evaluated at n = 2 since they become singular there. More details concerning this singularity 
will be discussed in the next section. 

The coefficients and V'^g in Table II take extremely large values at n = 6. The values 
of VT^ and V^f^ in Table III are also large. This is due to the fact that X>!! and have 
terms with the factors jz^f ^^^^ -^_]n , respectively, and that and d^g happen to be very 
close to one at n = 6. Actually we obtain dT^ = 0.995846 and = 1.00035 forn/ = 3 
(Table 11), and d'^^ = 1.07427 and d%=s^ = 1.08038 for n/ =4 (Table 111). But wc sec from 
(2.29) that and T>pfg are multiplied, respectively, by the factors 1 — ( f^psj j and 

I- y s —1-1 

^ ~ [ a (p2) j which become very small when and are close to one. Thus the 

contributions of the parts with 1)^=6 and T>'^^ to the 6-th moment of F^{x, Q^, -P^) do not 
stand out from the others. 



IV. SUM RULE OF F^{x, P^) 

The sum rule of the structure function F^, 

rdxF^{x,Q^P') , (4.1) 

^0 

can be studied by taking the n — 2 limit of Eq.(2.29). At n = 2 one of the eigenvalues of 
ln=2{g), the anomalous dimension matrix in the hadronic sector given in (2.13), vanishes, due 
to the conservation of energy momentum tensor. Thus we have a zero eigenvalue, — 0, 
for the one-loop anomalous dimension matrix ^^=2 ^nd, therefore, we get dT"^ — = 0. 
Among the coefficients which appeared in (2.29), two of them, namely, and E'^ would 
develop singularities at n = 2, since those coefficients have terms with the factor However, 

as we see from (2.29), both and are multiplied by a factor 
also vanishes at n = 2. Provided that we regard the expression -{1 — x^) as its hmiting value 



which 
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TABLE II: Numerical values of , ylf , fif , , , J:f (i = +,-,A^S'), and C" and Q'' for n = 
2,4, •• • ,12 in the case of Uf = 3. The calculation of and 2?jv5 performed by using 

the "exact" values of K^^g"", K^^'"" and K^^'"' given in Ref.[10] (at upper levels) and also by using 
the approximate expressions, i^^'^approx' -f^V'^approx and K^g^^to^ defined in Eqs.(3.24)-(3.26) (in 
parentheses at lower levels). 



n 


/•n 

^+ 




rn 


A1 




An 






"ns 


C" 


2 


0.4690 


0.4267 


0.4248 


-2.8403 




-5.5940 


1.7481 


-1.8535 


0.8290 


-9.3333 


4 


0.004336 


0.3639 


0.1836 


-0.5543 


-2.6267 


-1.3299 


0.07353 


3.3149 


1.4607 


-10.7467 


6 


0.0005428 


0.2324 


0.1164 


0.06133 


-1.8806 


-0.9403 


0.01652 


2.9783 


1.5349 


-9.1088 


8 


0.0001493 


0.1689 


0.08451 


0.009544 


-1.6566 


-0.8277 


0.006245 


2.9612 


1.4906 


-7.7504 


10 


0.00005803 


0.1318 


0.06591 


0.002817 


-1.5336 


-0.7664 


0.002993 


2.8263 


1.4169 


-6.7116 


12 


0.00002748 


0.1075 


0.05375 


0.001087 


-1.4425 


-0.7210 


0.001652 


2.6744 


1.3390 


-5.9074 



73" 



cn 




cn 




J7n 


'pn 


gn 


10.5867 




-10.9168 


6.9729 


-13.7973 


3.7817 


-251.3619 


-9.3990 


-23.9288 


-10.5807 


1.3106 


48.8620 


20.6599 


-204.5836 


1.8667 


-24.0991 


-12.4017 


0.4596 


56.4575 


29.3881 


-176.9466 


0.3993 


-29.0367 


-14.5993 


0.2217 


67.5380 


34.0579 


-157.4181 


0.1453 


-32.8877 


-16.4753 


0.1249 


73.0111 


36.6197 


-142.6108 


0.06532 


-35.8891 


-17.9598 


0.07764 


75.9024 


38.0024 


-130.8717 



60.5098 
( 60.5014 ) 

7.9871 
{ 7.9873 ) 

0.01877 
( 0.01877 ) 

0.03222 
( 0.03221 ) 
-0.001732 
{ -0.001738) 

-0.01825 
( -0.01826 ) 



32.9286 
( 32.9251) 
25.9791 
( 25.9222) 
-4007.0415 
(-4011.3304) 
165.7976 
165.9277) 
109.3285 
109.4090) 
86.8381 
86.9019) 



63.1965 
63.1909) 
11.6147 
11.5840) 
24025.6303 
24050.8845) 
82.2116 
82.2758) 
54.5447 
54.5847) 
43.3780 
43.4098) 



for e — > 0, —Inx , then the A1 and £^ parts of (2.29) give finite contribution as 



lim Al 

n^2 



lim 

n->2 



1 - 



1 - 



as{P^)) 
as{P^)) 



-r='in 



«.(P2) ' 



(4.2) 
(4.3) 
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TABLE III: Numerical values of -Bf, Pf, , (i = +,-,NS), and and a" for n = 

2,4, •• • ,12 in the case of nj = 4. The calculation of and 2?jv5 performed by using 

the "exact" values of K^^g"', K^^'"' and K^^'^ given in Ref.[10] (at upper levels) and also by using 
the approximate expressions, i^S'^approx' -f^V'^approx 
parentheses at lower levels). 



and K^G^^pro^ defined in Eqs. (3.24)- (3.26) (in 



n 


rn 




'-NS 


Al 


A1 


An 






on 

"ns 


C" 


2 


0.8078 


1.0582 


0.6231 


2.7608 




-6.0944 


3.8774 


-8.5894 


1.3076 


-16.3237 


4 


0.009356 


0.7327 


0.2661 


5.1244 


-3.7321 


-1.3858 


0.1688 


0.4820 


2.1599 


-18.7956 


6 


0.001235 


0.4656 


0.1679 


0.09529 


-2.9038 


-1.0480 


0.03909 


6.0485 


2.2.395 


-15.9311 


8 


0.0003465 


0.3374 


0.1215 


0.01953 


-2.7046 


-0.9735 


0.01495 


5.9573 


2.1613 


-13.5552 


10 


0.0001362 


0.2627 


0.09461 


0.006354 


-2.5904 


-0.9322 


0.007216 


5.6671 


2.0468 


-11.7384 


12 


0.00006497 


0.2140 


0.07704 


0.002598 


-2.4906 


-0.8963 


0.003999 


5.3501 


1.9293 


-10.3319 



XI" 



en 


£n 


pn 
'-'NS 








gn 


13.2519 




-12.7900 


7.0067 


-68.4928 


0.8275 


-439.6247 


92.4633 


-2.4550 


-11.2489 
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-28.0786 


24.5807 


-357.8108 


3.0154 


-37.7241 


-13.9828 


1.0159 


99.2476 


37.1197 


-309.4744 


0.8428 


-47.7463 


-17.3117 


0.5046 


121.0094 


43.9510 


-275.3197 


0.3366 


-55.8735 


-20.1683 


0.2888 


132.4677 


47.7946 


-249.4221 


0.1600 


-62.2711 


-22.4456 


0.1813 


139.2236 


49.9533 


-228.8909 



-84.4549 
(-84.4748 ) 
-17.3048 
(-17.3044 ) 
0.4163 
0.4163 ) 
-0.04747 
-0.04749) 
-0.2306 
-0.2306 ) 
-0.7548 
-0.7548 ) 



64.6182 
( 64.6102) 
-140.7574 
(-140.9078) 

894.6070 
( 895.0955) 

326.5791 
( 326.7480) 

228.9242 
( 229.0460) 

183.5002 
( 183.6024) 



63.5804 
( 63.5722) 
-64.7231 
( -64.7847) 
301.7867 
( 301.9492) 
116.7791 
( 116.8392) 
82.2373 
( 82.2809) 
65.9952 
( 66.0319) 
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The coefficient functions, anomalous dimensions and photon matrix elements at n = 2 
are given in Appendix D. Using these values we obtain ^4!! ^ = —1.3274 (—2.2857) and 
TT'^ = 5.7664 (18.553) for n/ = 3 (4). The numerical values of 

^n=2^ _^n=2 ^ + _ ^v^) ^ud ^^=2 and ^"=^ except for and E'T^, were already 
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given in Table II (for rif — S) and III (for Uf — A). 

Let us express the sum rule in the following form as 

dxF^ix, Q^ P') = ^:^{^^c^o + c^LO + ^^c^nlo + 0{al)] , (4.6) 

where the first, second and third terms in the curly brackets correspond to the LO, NLO 
and NNLO contributions, respectively. The coefficients clo, cnlo and cnnlo depend on 
the number of the active quark fiavors n/, and also on as{Q'^) and as{P'^). For the QCD 
running coupling constant q;s(Q^), we use the following formula which takes into account 

the j3 function parameters up to the three-loop level [41], 



47r (5^L {poLyPo (PoL)^ \Po 



2 



2 J Pi 4 



O 



1 

(4.7) 

where L^\n{Q'^ /A'^), and Po, Pi and P2 are given in Eqs.(3.2)-(3.4). Taking A = 0.2 GeV, we 
get, for example, a, (Q^ = lOOGeV^) = 0.1461 (0.1595) and a, (Q^^SGeV^) =0.2487 (0.2717) 
for the case = 3 (4) . 

We list in Table IV the numerical values of the coefficients clo, cnlo and cnnlo for the 
cases n/ = 3 and 4. We have studied three cases: ((3^ P^) =(30GeV^ IGeV^), (lOOGeV^ 
IGeV^) and (lOOGeV^ 3GeV2). We already know that Cnlo takes negative values [17]. We 
find that the coefficient cnnlo also takes negative values which are rather large in magnitude 
compared with those of c^o and cjvlo- Also listed in Table IV are the NLO {aag) and NNLO 
{aal) corrections relative to LO (a) and the ratios of the NNLO to the sum of the LO and 
NLO contributions for the sum rule of F2 {x , Q'^ , P^) . We see that the NNLO corrections 
give negative contribution to the sum rule. In fact, we will see in the next section that 
the NNLO corrections reduce F2{x,Q'^,P'^) at larger x. For the kinematical region of 
and which we have studied, the NNLO corrections are found to be rather large. When 
p2 = lGeV^ and Q^^SO^ lOOGeV^ or p2^3GeV^ and Q2 = 100GeV^ and Uf is three or 
four, the NNLO corrections are 7% ~ 10% of the sum of the LO and NLO contributions. 
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TABLE IV: The numerical values of coefficients clo^ cnlo and cmnlo in Ec[-(4-6), and the NLO 
and NNLO corrections relative to LO for the sum rule of F2 {xjQ"^ , P'^) in several cases of and 
P^. The ratios of the NNLO to the sum of the LO and NLO contributions are also listed. For the 
QCD running coupling constant Ug, we have used the formula given in Eq.(4.7) with A = 0.2 GeV. 





Q2(GeV2) P^{GeV^) 


CLO 


CNLO 


cnnlo 


LO 


NLO 


NNLO 


NNLO/(LO-FNLO) 


n/ = 3 


30 


1 


0.7631 


-11.66 


-331.2 


1 


-0.2063 


-0.0791 


-0.0997 




100 


1 


0.8613 


-12.21 


-355.3 


1 


-0.1649 


-0.0558 


-0.0668 




100 


3 


0.6690 


-11.22 


-313.8 


1 


-0.1949 


-0.0634 


-0.0787 


n/ = 4 


30 


1 


1.429 


-18.90 


-525.7 


1 


-0.1950 


-0.0800 


-0.0993 




100 


1 


1.614 


-19.59 


-551.4 


1 


-0.1541 


-0.0551 


-0.0651 




100 


3 


1.257 


-18.38 


-507.5 


1 


-0.1855 


-0.0650 


-0.0798 



V. NUMERICAL ANALYSIS OF F^{x,Q'^,P^) 

We now perform the inverse Mellin transform of (2.29) to obtain F2 as a function of x. 
The n-th moment is denoted as 

M^(n, Q\ = dx . (5.1) 

Jo ^ 

Then by inverting the moments (5.1) we get 

F^{x,Q^,P^) _ 1 i-^^'^ 



1 pLy+tOO 

— dnx--M^{n,Q',P') , (5.2) 

27^^ Jc-ioo 



X 

where the integration contour runs to the right of all singularities of (n, Q^, -P^) in the 
complex n-plane. In order to have better convergence of the numerical integration, we change 
the contour from the vertical line connecting C — ioo with C + ioo (C is an appropriate 
positive constant), introducing a small positive constant e, to 

n = C — e\y\ + iy, —00 < y < 00 . (5.3) 

Hence we have 

- / rRe{M2^(^,Q^p2)g-.in(x)-^ 

- elm{M]iz, Q\ P^)e-^i°(-)}] dy , (5.4) 
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FIG. 3: Virtual photon structure function F^{x,Q'^,P^) in units of (3an/(e^)/7r) ln(Q2/p2-) f^j. 
Q2 = 30 GeV^ and P'^ = l GeV^ with 7t,/ = 4 and the QCD scale parameter A = 0.2 GeV. We plot the 
Box (tree) diagram contribution including non-leading corrections (short-dashed line), the QCD 
leading order (LO) (dash-dotted line), the next-to-leading order (NLO) (long-dashed line) and the 
next-to-next-to-leading order (NNLO) (solid line) results. 

where z = C — ey + iy. 

As wc sec from Eqs.(2.29)-(2.37), the n-th moment M2(n,Q^,P^) is written in terms of 
coefficient functions, anomalous dimensions and photon matrix elements, which in turn are 
expressed by the rational functions of integer n and also by the various harmonic sums [42] . 
Thus we need to make an analytic continuation of these harmonic sums from integer n 
to complex n. There are several proposals for this continuation [43, 44]. The method we 
adopted here is to use the asymptotic expansions of the harmonic sums and their translation 
relations. The details are explained in Appendix E. 

In Fig. 3 we plot the virtual photon structure function F2{x, Q^, P^) predicted by pQCD 
for the case of n/ = 4, (5^ = 30 GeV^ and P^ = l GeV^ with the QCD scale parameter A = 0.2 
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GeV. The vertical axis corresponds to 



F^,{x,Q\P')l—nf{e')hi 



(5.5) 



p2 ■ 



Here we show four curves; the LO, NLO and NNLO QCD results and the Box (tree) diagram 
contribution including non-leading corrections. The box contribution is expressed by [17] 



where power corrections P^/Q^ and quark mass effects are ignored. It is noted that, in these 
analyses, even for the LO and NLO QCD curves, we have used the QCD running couphng 
constant as{Q'^) which is valid up to the three-loop level and is governed by the formula 
(4.7) and we have put A= 0.2GeV. 

The LO and NLO QCD results with the same values of n/, and as well as the Box 
contribution were already given in Fig. 6 of Ref.[17]. But in Ref.[17] the formula for asiQ"^) 
which is valid in the one-loop level was used to obtain the LO curve, while the two-loop level 
formula for asiQ"^) was applied for the NLO graph, and the QCD scale parameter A was set 
to be 0.1 GeV in both cases. The LO result in Fig.3 has a similar shape as the corresponding 
one in Ref.[17] but is different in magnitude; the former is shghtly larger than the latter for 
almost the whole x region. This is due to the fact that the one-loop-level formula for as{Q'^) 
was used for the LO curve in Ref. [17] while we applied the three- loop-level formula even for 
the LO result. On the other hand, the NLO curve in Fig.3 is similar to the corresponding 
one in Ref. [17] in shape and magnitude. 

Now we observe in Fig.3 that there exist notable NNLO QCD corrections at larger x. 
The corrections are negative and the NNLO curve comes below the NLO one in the region 
0.3 < X < 1. This is expected from the n — 2 moment analysis in Sec.IV. From Table IV we 
see that the ratio of the NNLO to the sum of the LO and NLO contributions for the sum 
rule of F:^{x, P^) is -0.099 for the case of Uf = 4, = 39 GeV^ and P^ = 1 GeV^. At 
lower X region, 0.05 < a; < 0.3, the NNLO corrections to the NLO results are found to be 
negligibly small. 

We have also studied the QCD corrections to P^(x, Q^, P^) with different and P^ but 
with n/ = 4. In Fig.4 we plot the case for = 100 GeV^ and P^ = 1 GeV^. Another case 
for (5^ = 100 GeV^ and P^ = 3 GeV^ is shown in Fig. 5. We have not seen any sizable change 




(5.6) 
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for the normalized structure function (5.5) for these different values of and P^. In both 
cases the NNLO corrections reduce Q^, -P^) at larger x. We have examined the n/ = 3 

case as well. It is observed that the normalized structure function (5.5) is insensitive to the 
number of active flavors. 

Finally we should note that the spin-averaged structure function directly accessible in 
the experiment is not F^^ but the so-called effective structure function ~ F2'^ -|- (3/2) F2 
as discussed in ref. [15, 21, 45] 



VI. LONGITUDINAL STRUCTURE FUNCTION F2(x, Q'^ , P^) 

We have considered the structure function so far. Regarding another structure 
function F2, its LO contribution, which is of order a, was calculated in QCD for the 
real photon (P^ = 0) target in Refs.[5, 6]. The analysis was extended to the case 
of the virtual photon (A^ -C P^ -C Q^) target [17]. We will now derive a for- 
mula for the moment sum rule of F2{x, Q^, P^) up to the NLO {0{aas)) corrections. 
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X 

FIG. 5: Virtual photon structure function F^{x,Q^,P^) for Q^^iqo GeV^ and P^ = 3 GeV^ with 
n/ = 4 and A = 0.2 GeV. 



Comparing Eq.(2.8b) with (2.8a) and examining the form of (2.11), we see that the 
formula for F2{x,Q^, P'^) is obtained from (2.18) only by replacing C2,n(l, ^(Q^)) and 
C2_„(l, ^(Q^), a) with Ci^„(l, ^(Q^)) and C2 „(1, ^(Q^), «), respectively. An expansion is 
made for Cl^„(1, ^(Q^)) and CI ^{l,g{Q'^), a) up to the two-loop level as 

CMQ'^)) = cZ + ^c»i + ggc£ + . . . , (6.1) 

S(Q%a) = ^C2« + + . . . . (6.2) 

Here we note that there is no contribution of the tree diagrams to the longitudinal coefficient 
functions and thus we have C^l \ — 0. 

The moments of Fl{x,Q'^ , P"^) are then given as follows (see Eqs.(2.29)-(2.37) for com- 
parison) : 
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(6.5) 



pn;^(l) pn 



In n / J \n _ \n 



A" - A" + 2/3o 
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(6.6) 



1 



pn- (1) pn 

2/3o(cz!?+i:«-ca), 



(6.7) 
(6.8) 



with i,j = +,—,NS. The coefficients 'B"^.) i ^'^'^ ^(1) represent the LO terms [5, 6, 17], 
while the terms with S^^^-, J'J^m and are the NLO {aag) corrections and new. It is 
noted that, among these coefficients, £'^^^ _ becomes singular at n = 2 since it has terms 
with the factor ^ and vanishes as n— >2. But again, as in the case of the moments of 

r / \dri 

F2 {x , Q"^ , P'^) , this coefficient is multiplied by a factor 1 — ias{Q^)/ois{P^)] , and thus 
the product remain finite at n = 2. 

The one- loop longitudinal coefficient functions were well known [28, 46, 47]. They are 
written in the form as 



.G(l) 



NS{1) 



■y7(1) 
L,n 



(6.9) 



where S^^^ = Bns,l^ Pg,l ^^'^ -^"l given, for example, in Eqs.(6.2) - (6.4) of Ref.[6]. 
The two-loop longitudinal coefficient functions corresponding to the hadronic operators were 
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calculated in the MS scheme in Ref.[31, 32] ^. The results in MeUin space as functions of n 
are found, for example, in Ref. [33]: 

ct^ = {c^^ri^) + c^i!ri^)}, (6.10) 

= c?i(n), (6.11) 

Cl'J'' = S.s cf-\n) . (6.12) 

where ci^^^{n), ci^^^{n) and c^l\{n) are given in Eqs.(203), (204), and (205) in Appendix 
B of Ref. [33], respectively, with N being replaced by n. The two-loop photon longitudinal 
coefficient function C^^^ is expressed as 

C'S = ^T^ai^) , (6-13) 
and c^l\{n) is obtained from c^l\{n) in (6.11) by replacing — > and ^ — > 1. 

Inverting the moments (6.3), we plot in Fig.6 the longitudinal virtual photon struc- 
ture function F2{x, Q^, P^) predicted by pQCD for the case of n/ = 4, = 30 GeV^ and 

= 1 GeV^ with the QCD scale parameter A = 0.2 GeV. The vertical axis is in units of 
F2{x, Q^, P^)/^n/(e^). Here we show three curves; the LO and NLO QCD results and the 
Box (tree) diagram contribution, which is expressed by 

F2'^''°''\x,Q\P^) = ^nf{e^)[4x\l-x)^ . (6.14) 

The LO result in Fig.6 is consistent with the corresponding one in Fig. 5 of Ref. [17], al- 
though the formulae used for as{Q'^) differ in detail. We see from Fig.6 that the NLO 
QCD corrections are negative and the NLO curve comes below the LO one in the region 
0.2 < a; < 1. 

The QCD corrections to F2{x,Q^ , P^) for different values of Q^, and Uf are also 
studied. The case for Q"^ = 100 GeV^ and P^ = 1 GeV^ with Uf = A is shown in Fig. 7. 
The LO curve has hardly changed from the one for Q"^ — 30 GeV^ and P^ = 1 GeV^. The 
NLO corrections get smaller. The LO and NLO QCD curves for Q^^lOO GeV^ and p2 = 3 
GeV^ with Uf — A appear to be almost the same as those in the case of (5^ = 30 GeV^ and 
P^ = l. The cases for nf = 3 are examined as well and we find that the normalized function 
F2{x, Q'^, P^)/^nf{e'^) is insensitive to the number of active flavors. 

^ The earlier calculations [48-51] were found partly incorrect. For quark coefficient functions c^^^^^^ and 
^L^q^"" Eq.(6.10), there is a complete agreement between Ref. [51] and [31, 32, 52] while for gluon 
coefficient c'l^ in Eq.(6.11) the result of Ref.[49] was corrected in Ref. [53]. 
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FIG. 6: Longitudinal photon structure function F2{x,Q^,P'^) in units of {3an f {e'^) / tt) for Q^ = 30 
GeV^ and P'^ = l GeV^ with n/ = 4 and the QCD scale parameter A = 0.2 GeV. We plot the Box 

(tree) diagram contribution (short-dashed line), the QCD leading order (LO) (dash-dotted line) 
and the next-to-leading order (NLO) (long-dashed line) results. 

VII. CONCLUSIONS 

We have investigated the unpolarized virtual photon structure functions F^{x,Q'^, P'^) 
and F2{x, Q^, P^) for the kinematical region <C <^ Q'^ in QCD. In the framework of 
the OPE supplemented by the RG method, we gave the definite predictions for the moments 
of P^(a;, Q"^, P^) up to the NNLO (the order aas) and for the moments of F2{x, Q'^, P^) up to 
the NLO (the order aag) . In the course of our evaluation, we utilized the recently calculated 
results of the three-loop anomalous dimensions for the quark and gluon operators. Also we 
derived the photon matrix elements of hadronic operators up to the two-loop level. 

The sum rule of F^ix, Q^, P^), i.e., the second moment, was numerically examined. The 
NNLO corrections arc found to be 7% ~ 10% of the sum of the LO and NLO contributions, 
when P2 = IGeV^ and Q^ = 30 lOOGeV^ or P^ = SGeV^ and = 100GeV^ and n/ is 
three or four. 
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FIG. 7: Longitudinal photon structure function F2{x, Q^,P^) for = iqo GeV^ and = 1 GeV^ 
with n/ = 4 and A = 0.2 GeV. 



The inverse Mellin transform of the moments was performed to express the structure 
functions F2{x,Q'^, P^) and F2{x,Q'^, P^) as functions of x. We found that there exist 
sizable NNLO contributions for at larger x. The corrections are negative and the NNLO 
curve comes below the NLO one in the region 0.3 < x < 1. At lower x region, 0.05 ^ x < 0.3, 
the NNLO corrections to the NLO results are found to be negligibly small. Concerning F2, 
the NLO corrections reduce the magnitude in the region 0.2 < a; < 1. 

The comparison of the present NNLO theoretical prediction for the virtual photon struc- 
ture functions with the existing experimental data will be discussed elsewhere. 
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APPENDIX A: EVALUATION OF Mn{Q^/P^,g{P'^)) AND Xn{Q^/P^,9{P^),a) 



In order to evaluate the integrals for Mn{Q'^/P'^,g{P'^)) given in (2.22), we employ the 
same method that was used by Bardeen and Buras in Ref.[6] and make a full use of the 
projection operators obtained from the one-loop anomalous dimension matrix j^: 
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where Xf {i — NS) are eigenvalues of 7^ and are expressed as 



"^NS — INS ' 



and Pp are the corresponding projection operators, 
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With an expansion of P{g) up to the three- loop level in (2.24), we get its inverse as follows: 
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Then using (Al) and (A6), we perform integration in (2.22). 
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The result is 
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(A7) 



where gi — g^{P^) and gi — g^{Q^)- The first term is the leading, and the second and third 
terms are the next-to-leading terms. The rest are the next-to-next-to-leading terms. 

Once we get the above expression for M„(QV^^^?(^^)) expanded up to the NNLO, 
we use an expansion of Kn{g,a) in (2.25) up to the three-loop level and can evaluate 
Xr,{QyP^,g{P^),a) in (2.17) up to the NNLO. The result is 
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where gf = g'^{P'^) and = g'^{Q'^). The first term is the leading, and the second to fourth 
terms are the next-to-leading terms. The rest are the next-to-next-to-leading terms. 
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APPENDIX B: S^sT°^(n), E^P'°''{n) AND Ep,^{n) 

We give the exphcit expressions of E^^^^°^{ri), E^^'^°^{n) and Epsjin) which have appeared 
in (3.24)- (3.26). They are obtained by taking the MeUin moments of the parameterizations 
for P^sl{x) and Pg^\x) and of the exact result for P^s\{x), which are presented in Eqs.(6)-(8) 
of Ref.[ll]. Using a single harmonic sum Sm{n) defined by 



128^1 (n)^ ^ 62.5244^1 (n)=^ 50.08^1 (n)^ 256S2{n)Si{nf 

27n n n + 1 9n 
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lbQ.2AS2{n)Si{n) _ imASri{n)Si{n) _ l2d>S2{nf ^ 785.14,gi(n) 

n + 1 27n 9n n 

325.45i(n) 300.485i(n) 175.3^2(^1) 520.8^2(^1) 150.2452(n) 

[n + 1)^ n v? {n + 1)^ 

591.15153(n) _ 100.16S3(n) _ 2565'4(n) _ 492.087 ^ 1262 _ 449.2 

n n + 1 9n n n + 1 n + 2 

1445 1279.86 1169 403.2 160 300.48 512 

n + ^ {n + iy~ ^ ^ ~ {n + 1)"^ ~ 9nP 

r 32Si{ny 258.U2Si{ny 270Si{ny 269.45i(n)2 

f \ — ~t~ ~t~ 

[ 27n n n + 1 v? 

535.244^2(^)^1(71) 905.065i(n) 5405i(n) 17.0465i(n) 




(Bl) 



they are expressed as 





n n+1 n^ n 



628.124 114.4 24.86 53.39 49.5895 26.63 




n n + 1 n + 2 n + 3 n^ {n + iy 
21.984 540 64 1 



(B2) 
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approx , 
'G7 I 



n 



rfa;a;"'^{thc r.h.s. ofEq.(7) in Rcf.fll]} 
2>2Si{nf 325i(n)3 79.13Si(n)2 79.13^i(n)2 433.25i(n)^ 



27n 27(n +1) n 
429.644^i(n)2 862.84452(n)5i(n) 



(n + l)2 
1512.395i(n) 



n 



+ 



1512.395i(n) ^ 549.55i(n) 



n + 1 



n + 1 
862.844^2 (n)5i(n) 
n + 1 

707.765i(n) 



2460^1 (n) 



4185.69,5i(n) ^ 628.63^2(ri) _ 628.63,S2(n) 2893.2,52(n) 



3756.0452H 



n 



n + 1 



(n + 1)^ n 
73.1409 1673.57 3180.43 



3324.03,53(71) ^ 3324.03,53(71) 



n — 1 
6400 



37l2 

358 



3688.39 
{n + 1)2 ~ 
/2460 



+ 

2247.4 



n + 1 
1420 406.7 



566.7 



- 1 
+ 



n + 2 
990.14 

(71+1)2 



n + 3 



+ 



128 



1600 



n + 4 ■ 3(71-1)2 
9438.76 3584 



371^ 



(71+1)^ 



97l5 



2460 \ 



9(71+1)5 V 71 

[32^1(71)2 32Si{ny 



/2460 



2460 



:)C2 



9.1335i(n) 



9n 



9(n + 1) 



18.26652(71)^1 (7z) 46.42645i(7i 



77. + 1 

23.29115i(7z) 



+ 



n 



(71+1)2, 

9.1335i(7z)2 

^ (77+1)2 

46.4264^1 (77) 

77. + 1 



18.266,S2(7i),5i(7i) 



+ 



77. 

16.185i(7z) 

n2 



76.665i(77) 113.192^1(77) 19.7356^2(77) 19.7356^2(77) 



(71+1)^ 



+ 



(71 + 1)2 



+ 



n 



n+1 



85.793^2 (ri) ^ 104.059,^2 (n) _ 94.926,^3(^1) ^ 94.926,^3(71) 



+ 



40.5597 

n — 1 
161.767 

(71+1)2 



(71+1)^ 



21.1683 17.0286 

h - 



77 
52.82 



77'^ 



+ 



77. + 1 

13.3489 



77 

93.37 

" 77 + 2 
128 



+ 



101.05 

77 + 3 

299 



77+1 

44.1 

77 + 4 



+ 



115.341 



77^ 



[n + 1)3 977^ (77 + ly 



(B3) 



Ep,^{n) = - [ rfxx"-^{the r.h.s. of Eq.(8) in Ref.[ll]} 
Jo 



. 2464 432 72 _ 38360 _ 344 _ 368 

~ ^ ' ~S1(77- 1) ^ ~ ^ 77 + 1 ~ 81(77 + 2) ~ ^ ~ (77. + 1) 

3584 288 208 448 96 96 

+ ^ + 7 7T^ + 



2 



27(77 + 1)2 773 (77 + 1)3 9(n + 2)3 774 (77 + 1)^ 
256 64 128 

'3(77 + 2)4 + ~ (77+1)5 



256 64 128 , 

+T77-^ + -^-7-^ • (B4) 



Note that Eps^{n) is an exact result. 
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APPENDIX C: MELLIN MOMENTS 45(n) AND WITH i = 1,2 

AND 4?(") AND bfj{n) 

The expressions of aqg{n) and bqg{n) with i —1 and 2 are obtained by taking the moments 
of the functions aqg {z) and bqg {z) as 



/" dzz''-^afg{z) , hfg{n)= f dzz''-^hfg{z) , with i = 1,2 (CI) 
JO Jo 

where aql{z) and hql{z) are extracted from the e-independent terms of A^^^^i^z,^,^ 
given in Eq.(A7) and of ^^^^{z,^,-^ in Eq.(A8) of Ref.[40], respectively. See also 
Eqs.(2.27) and (2.28) of Ref.[40]. 

The one-loop results are 

oSW - ?4|fi-^ + ^)fs,W-A.i,-^ + ^|,(C2) 



' 2 \\n n+l n + 2/\ ^ ' ) rv" (n + l)^ (n + 2)2 ^ 

= ^lef^ —\ . (C3) 

^ 2 Vn+1 n + 2/ ^ ^ 



"19 



The two-loop results are 

- ^-Y{(s-;r^+;r^)(-5^'W'-^^^(")^-(") + f^»(") 

-16>S2,i(n)-48C3) 
2/6 _ 24 32 _ 16 40 _ 32 \ 

/12 8 40 16 128 128 

\n ^ n+l ~ n + 2 ^ ^ ~ (n + 1)^ ^ (n + 2)2 

32 176 128 



7 + 


n + 2 


n-Fl 


8 


40 


n-Fl 


n + 2 


8 


16 


+ - 





3 ^n + l)3 (n + 2)3 



/6 8 16 16 40 32 \ 

+52(n) - + - — + ^^^-^^2 - („ + 2)2j 

14 38 48 64 126 80 22 

~ n + l ^ n + 2 ^ ^ ~ (n + l)2 ^ (n + 2)2 ~ ^ 
88 128 20 



(n + l)3 (n + 2)3 n^ (n + 1)^ 
n f / n f \ 2 
+ terms proportional to C^-^ V 2 / ' 
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11 12 10 



-n n+1 n + 2 (n + 1)^ (n + 2)^ 

3 4 1 1 9 8 6 

— + r -^ + ^ + 



n n + 1 n + 2 n^ {n + \f {n + 2f {n + Vf 



+ terms proportional to ^A-^ or i (C5) 

where S2,\{n) = Yl]=i ^'S'i(j). The terms proportional to (^"^j i^i cifaif^) bfgi.'"^) come 
from the external gluon self-energy corrections and should be discarded for the photon case. 



Similarly the expressions of agg{n) and hgg{n) are obtained by taking the moments 
of the functions afg{z) and bfg{z) which are extracted from the e-independent terms of 
A^g^^^l^z, given in Eq.(Al2) and of 7) in Eq.(Al3) of Ref.[40], respec- 

tively See also Eqs.(2.34) and (2.35) of Ref.[40] ^. 

The two-loop results for a^gj (n) and b^gj (n) are 



16 48 32 128 32 24 



+Si(n) — — + + — — — + ^ + 



9(n-l) n n+1 9(n + 2) (n + 1)^ 

64 32 48 



3(n + 2)2 n3 {n + lfJ 
/ 32 32 32 32 \ 

680 48 16 184 56 56 256 44 

+ 77^7 7T + 



3 



27(n-l) n n+1 27(n + 2) n^ {n + lf 9(n + 2)2 n 
76 128 40 88 55 1 

^(n + l)3 " 3(n + 2)3 ~ ^ ' {n + ly ' Y ^ ^'j 

-I- terms proportional to or (~^) ' ^^^^ 

Two terms Jggbl'g^^ and Jgqbqg^^ are missing in the e-independent terms of Eq.(2.35) of Ref.[40]. The 
both are needed in order to extract bgg {z) correctly. 
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,f2W ^ ^ I o / a/ 32 32 64 32 

b'-,J(n) = CF—<Si(n){— ^— - + - 

' 2\'^^V3n-l) n 3 n + 2 {n + If 



128 64 448 32 96 128 96 

+ 



9(n-l) n 9{n + 2) {n + lf 3(n + 2)2 {n+l] 



+ terms proportional to or C^-^ . (C7) 

The terms proportional to ^"^^ in Ogg^(n) again come from the external gluon self-energy 
corrections and should be discarded for the photon case. Furthermore, the contribution of 
the last two terms in the curly brackets of (C6), more exphcitely, Cp^ -'-^Cs^ > is also 

resulted from the external gluon self-energy corrections and is thus irrelevant. 
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APPENDIX D: VALUES AT n = 2 



1. Coefficient functions 



With 5^ — (e^), Sjsfs — i and 5^ — 3nf{e'^), we have 
• At tree-level 



^V(O) _ c ^G(O) _ p, ^NS{0) _ r ^7(0) _ r, 

^^2,71=2 — "V" ' ^^2,71=2 — ^ > <^2,n=2 " "iV5 , <^2,n=2 — ^ ■ 



At one-loop 



At two-loop 



2. Anomalous dimensions 



• At one-loop 

(0),n=2 _ (0),n=2_^ /16\ (0),n=2 _ / 4\ 

(0),n=2 ^ / 16 \ (0),n=2 /4\ 

and 
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At two-loop 



^,(l),n=2 

Ins 




128 
27 




— '-'A'-^F 27 / ^F^f \ 


208 
27 


(l),n=2 


= ^^^/( 27) +^^^/( 


148 
27 


(l),n=2 




■208 
. 27 


(l),n=2 

7gg 







224^ 
224^ 



(D6) 



and 



= 3n,(e^)C,(^) , (D7) 



• At three-loop 

We get from [26] and [27], 

6256 128 . \ o /41840 128 . \ _ ^ / 448 



,^2^/ 17056 i2o.\,^2 / 6824 128 X / 1120 256 x 

(2)n=2 /278 208 \ ^ 2/2116\ . / 3589 \ 

^ 2/568\ ^2 /14188 256 ^\ ^3 /1120 256 ^ \ 
+C,n}{-) + - -C3) + - -C3) . 

(2)n=2 / 278 208 ^\ ^2/ 2116\ ^2 /3589 \ 

= C^C^n/(-— + — Cs) + CAn}[-^) + C^>/(^ - 48(3) 

^ 2/346\ ^2 /4310 64 \ 
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and from [10] 



^(2),n=2 



8620 128 
" 243 ^ 

o /2vr^ / 692\ 6044 64^ \ ( 8620 128^ 

-1880\ ^ ^ / 1138 64 \ ^9 /9592 128 



(D9) 



Note that a relation 7^*^'" ^7gg" ^ ~ 7v;g" ^7Gvi" ^ = is indeed holds for i = 0, 1, 2. 



3. Photon matrix elements of quark and gluon operators 
• At one-loop 

(1)G 



At two-loop 



= , (DIO) 
I«f = 3n,((e'>-(.V)(?). 



= 3n,(e=)C,(|5) , (DID 
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APPENDIX E: ANALYTIC CONTINUATION OF THE HARMONIC SUMS 



The moments of F2{x, Q^, P^) given in (2.29) are expressed by the rational functions of 
integer n and the various harmonic sums. The single harmonic sums are defined by 

^»(")-E^. (El) 

where /c = ±1, ±2, • • • , and the higher harmonic sums are defined recursively as 

Sk,mi,-,mp{n') — Tj-^i 'S'^i,-,mp(j)) (E2) 

3=1 

where indices k and mi, • • • ,mp take nonzero integers. In order to invert the moments so 
that we get as a function of x, we need to make an analytic continuation of these harmonic 
sums from integer n to complex n. Since the moment sum rules of the s-u-crossing-even 
structure function are defined for even integer n, the continuation should be performed 
from even n. Thus whenever a factor (—1)" appears, it should be replaced by (+1). The 
method we adopted here for the analytic continuation is to use the asymptotic expansions 
of the harmonic sums and their translation relations. Choosing the following two harmonic 
sums 

n 

SM-Y.-,^ (E3) 



s 



i=i j=i ■> fc=i 1=1 



as examples, we explain how we get approximate analytic formulae for these sums. 
The asymptotic expansion of Si{n) for large n is well-known: 

The right-hand side has a simple analytic property. On the other hand, Si{n) satisfies the 
following translation relation: 

S,{n) = S^{n + l)--^ . (E6) 

This relation is valid not only for integer n, but also for complex n. Therefore, our algorithm 
to evaluate Si{n) at arbitrary complex n is as follows: (i) If \n\ > no, where no is some 
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positive integer at which the asymptotic expansion (E5) holds at a desired accuracy, then 
we use the expansion (E5) to evaluate Si{n). (ii) For \n\ < no, we apply the translation 
relation (E6) and shift the argument n ^ n + 1 repeatedly, until the shifted new n satisfies 
the condition \n\ > uq so that the asymptotic expansion (E5) for Si{n) may be used with a 
desired accuracy Then Si{n) is evaluated by a formula 

n—n ^ 

Si{n)^S,{?i)-J2-—:, (E7) 

1=1 

where the expansion (E5) is used for Si{n). 

In the case of a more complicated higher harmonic sum 5'i,i,_2,i(?T') with even integer n, 
its asymptotic expansion for large n is given by 

SlT-2 lin) = Co,.ln^(n) + Co,ln(n) + Co,o + '-^^ + ^ + '-^^ + ^ + • • • , (E8) 
' ' ' n n n-^ n-^ 

where 

co,2 = ci,i = -^C(3) = -0.37564278- • • , 

co,i = -^7^C(3) - ^C'(2) = -0.63658940 • • • , 

Co,o = -|^7^C'(2) - ^7lC(3) - ln(2)Li4 - ^ln^(2)C(3) + ^ln3(2)C(2) 

-lln5(2) + ic(2)C(3) + ^C(5) - Lis Q) = -0.89930722 • • • , (E9) 

ci,o = -^7i.C(3) - ^C'(2) + ^C(3) = 0.057348080 • • • , 

C2,i = ^C(3) = 0.062607130- ••, 

= —7i5C(3) + —C^(2)-—C(3) = -0.22868296- •• . 
96 ' ^ ' 160 ^ ' ' 

Also 'S'^™"_2,i(^) satisfies the following translation relation: 

ST'^'^ oAn)^ ,5rr 2 1 + 2) - I + iin + 2) + ^ -S^'Ti in + 2) . 

i,i,-2,iv J i,i,-2,iv ; \^n+l n + 2j ^'"^'^^ ' (n + l)(n + 2) "^'^^ ' 

(ElO) 

Note that the right-hand side of (ElO) is written in terms of the same harmonic sum S'J^'ig i 
and lower harmonic sums S'J™2 1 3^2^ but with a larger argument n+2. When \n\ > Uq, 
the asymptotic expansion (E8) is used to evaluate 'S'i^f'i2,i(^)- For \n\ < no, we apply (ElO) 
and shift the argument n ^ n + 2 repeatedly, until the shifted new n satisfies the condition 
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\n\ > no so that the asymptotic expansion (E8) of Siy^_2^i{n) may be used with a desired 
accuracy. The lower harmonic sums 'S'l^f 2,i(^) ^^i^i''^) evaluated in a similar fashion. 

In practice, we take Uq = 16 for all harmonic sums. Then the asymptotic expansion 
formula for each harmonic sum is derived so as to ensure double precision accuracy (15 
significant figures) at n = no. For example, Si{n) and 'S'^^^^_^2,i(^) expanded up to the 
terms with l/n^° and respectively. 
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